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Abstract 

We investigate the structure and mobility of dislocations in hep 4 He crystals. In addition to fully 
characterizing the five elastic constants of this system, we obtain direct insight into dislocation 
core structures on the basal plane, which demonstrates a tendency toward dissociation into partial 
dislocations. Moreover, our results suggest that intrinsic lattice resistance is an essential factor in 
the mobility of these dislocations. This insight sheds new light on the possible correlation between 
dislocation mobility and the observed macroscopic behavior of crystalline 4 He. 



1 



Dislocations are line defects that play a central role in the mechanical deformation be- 
havior of crystalline solids [1]. Their activity is widely known to be pivotal in classical 
solids, controlling phenomena such as fracture and brittle-ductile transitions in metals and 
semiconductors. Much less is known [2], however, about the influence of dislocations on 
the properties of quantum crystals, which are solids where the quantum-mechanical zero- 
point kinetic energy is significant compared to the typical energy scale of the interatomic 
interactions JjJ. 

Recent experimental studies on crystalline solid 4 He, the propotypical quantum solid, 
indirectly indicate that dislocations are indeed involved in macroscopic phenomenona. The 
occurrence of apparent superfluidity j^, 5| and the observation of elastic stiffening |s 8 1 , 
for instance, have been linked to the mobility of dislocations. However, in contrast to 
classical solids, for which an abundant body of experimental results exists, the lack of specific 
experimental data concerning the behavior of dislocations prevents a direct investigation of 
this relationship. In view of these difficulties, one needs to resort to theoretical approaches. 
In principle, a realistic picture of dislocations and their properties is possible using path- 
integral or variational Monte Carlo simulations. However, given that dislocation modeling 
requires a simultaneous treatment of different length scales [9j, one associated with the core 
properties and the other associated with distances large compared to the atomic scale, the 
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viability of these direct approaches remains a subject of debate 

In this Letter, adopting a methodology that incorporates both of the relevant scales, we 
obtain fundamental insight into the properties of dislocations in the archetypal quantum 
crystal: solid hep 4 He. We compute the intrinsic dislocation structure and mobility of four 
different dislocation types using the multiscale paradigm of the semi-discrete Peierls-Nabarro 
(PN) model 13l-|l5|. It constitutes a hybrid continuum-atomistic approach that captures 
the long-range elastic fields as well as the lattice-discreteness effects associated with the 
dislocation core. All parameters in the model are determined using quantum-mechanical 



expectation values for 4 He applying the shadow wave function (SWF) formalism 16|, Il7| . 
In addition to providing key information concerning the elastic properties, our results shed 
new light onto the possible role of dislocations in the experimental observation of the elastic 
stiffening of solid and its potential connection to the apparent superfluidity in this quantum 
crystal. 

Within the semi-discrete PN framework, a straight dislocation lying along the y-direction 
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is represented in terms of a set of misfit vectors 5, that describe the disregistry of atomic 
row i relative to their counterparts on the other side of the glide plane, as depicted in Fig{T] 
Panel a) shows a schematic representation of two sets of atomic rows extending along the 
y-direction, one set on each side of the glide plane. The rows above the plane are labeled 
by the index i. Panel b) shows a top view, depicting a disregistry vector c\ for row %. In 
this manner, the total misfit associated with a given dislocation, described by the Burgers 
vector b, is thought of as distributed among the atomic rows along the x-axis, subject to the 
boundary conditions <5_oo = and 5^ = b . Here, we consider two-dimensional misfit vectors 
5i = 5fx + 5fy, with edge and screw components along the x and y directions, respectively 
The equilibrium structure of the dislocation is then represented by that particular misfit 
distribution Si that minimizes the dislocation energy per unit length, 

Udisl ^elastic ^misfit ^stress + c, (1) 

where 



Elastic = Y^Xij (K e P\ P) + K s P l P fj , (2) 

i,3 

t/misfit = y^j{5i)Ax, (3) 

i 

Stress = hJ2 i re & + 5 ti) + t s (5! + SU)] Aa:, (4) 



and C is a constant that can be ignored 



13] . In the elastic part of the dislocation energy, 



Eq. ([2~lh Xij is a discretized universal kernel [14J |. and K e = fi/4n(l — v) and K s = fi/Air 
are elastic pre-factors with \i the shear modulus and v the Poisson's ratio. In addition, 
p\ = (5f — Sf_j) /As and pf = (5 4 s — /Ax, where Ax is the distance between adjacent 
rows in the defect-free crystal. Eq.([3]) represents the misfit contribution, in which ^(5) is 



It describes the excess 



known as the generalized stacking-fault (GSF) energy surface 
energy per unit area of a crystal that is subjected to the following procedure. It is first cut 
into two defect-free parts across a given plane. The two parts are then displaced relative 
to each other by a vector 5, after which they are patched together again. An example 
configuration of the GSF on the basal plane of the hep structure is the intrinsic stacking fault 
(ISF), in which the displacement vector 5 describes the associated shift in the planar stacking. 
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In the context of the PN model, the GSF surface reflects the inter-atomic interactions in the 
system and serves to model the details of the dislocation core on the atomic scale. Finally, 
the stress term of Eq. (jlj) accounts for the work done by any external stresses, where r e and 
r s denote the magnitude of the components of the stress tensor that couple to the edge and 
screw displacements, respectively [l]. The quantities that specify the model for dislocations 
in a particular material are the elastic parameters /i and v , the GSF surface 7(<5) associated 
with the glide plane of interest, and Ax. 

Here, we employ the SWF model based on the parameter set of Ref.{l3| to determine these 
quantities for solid hep 4 He (space group 194). In order to determine its elastic properties, 
we employ a computational cell containing 720 particles at a density of 0.0294 A~ 3 , which 
corresponds to lattice parameters a=3. 63668 A and c=5. 93866 A, subject to standard pe- 
riodic boundary conditions. Sampling configurations according to the quantum-mechanical 
probability density of the SWF model using the Metropolis algorithm, we then compute 
expectation values of the stress tensor [19] associated with the six independent deformations 
of the periodic cell, imposing strain levels of 0.25%. Using the standard relationship between 
the stress and strain tensors 2(|, we extract the five independent elastic constants of the 
hep structure. The results, which, to the best of our knowledge represent the first complete 
estimate of the elastic constants in hep 4 He, are reported in Table [H The shear modulus 
/i = C44 = 17.1 ± 0.8 MPa, is in good agreement with the value of 14 MPa that follows from 
the ratio of the experimental shear stress and strain values reported in 2l|. Poisson's ratio, 
obtained from the results in Table (TJ is found to be v = 0.151. Both theoretical values are 
those corresponding to shear directions in the basal plane, in which hexagonal crystals are 
isotropic 22 1. 

Since basal slip is known to be the dominant dislocation glide mechanism in hep solid 4 He 



[23] we focus on the properties of these particular dislocations and compute the GSF surface 
associated with the basal plane. For this purpose, we utilize the 720-atom cell and impose a 
series of 400 slip vectors 5 in the basal plane by adjusting the periodic boundary condition 
along the c-axis. The shadow degrees of freedom of the atoms immediately adjacent to the 
slip plane are allowed to vary only along the c-direction to maintain the relative displacement. 
Using the Metropolis algorithm we then sample configurations according to the SWF and 
compute the expectation value of the Hamiltonian as a function of 5. Subtracting the 
expectation value at 6 — and dividing by the area, we then obtain the GSF surface. 
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In order to implement the results into the PN mode, we fit the results using a Fourier 
series that reflects the lattice symmetry of the basal plane of the hep structure: 7(5) = 
^]c ( jexp(iG • 5), in which we use a set of 81 two-dimensional reciprocal lattice vectors G. 

G 

The results are shown in Fig. [2] The perfect crystal configuration, which has GSF value 
of zero, is associated with the displacement 5 = and its periodic equivalents. The ISF 
configuration, which corresponds to the displacement vector (and equivalents) 5 = (0, b p ), 
with bp = ^\/3a = 2.0996 A the length of a partial Burgers vector, has an excess energy of 
0.0063 mJ/m 2 . 

Using our estimates for the elastic properties and the GSF surface in the PN model, we 
investigate the structure and intrinsic mobility of 4 dislocation types on the basal plane: (i) 
screw, (ii) 30°, (iii) 60°, and (iv) edge. Fig. [3] shows the optimized disregistry profile for the 
screw dislocation, obtained by minimizing Eq. ([T|) at zero external stress. As expected, given 
the low stacking-fault energy (SFE) value , it dissociates into two 30°partial dislocations with 
opposite edge components, separated by an ISF area with a width of 29 atomic rows, which 
corresponds to 91.33 A. The core width <j of the partials, defined as the distance over which 
the displacement changes from | to | of it total value 14j, is approximately 1 atomic row 
or ~3.1 A. The second line of Table HT1 contains the dissociation widths of the other three 
dislocations, showing an increasing ISF width with increasing edge component, consistent 
with dislocation theory [l|. 

In addition to the structural properties described here, the PN model permits an estimate 
of the intrinsic dislocation mobility, which is measured in terms of the Peierls stress. To 
this end, we impose an external shear stress in the glide plane parallel to the total Burgers 
vector. It produces maximal force per unit length [lj on the dislocation line for the given 
stress magnitude. This magnitude is then increased in small steps, followed by minimization 
of Eq. (CQ) with respect to the disregistry vectors 5{ . At a critical stress value, the so-called 
Peierls stress, an instability is reached and an equilibrium solution ceases to exist. 14, Ijj] 
In this situation the dislocation becomes free to move through the crystal. The third line 
of Table [III contains the Peierls stress values for the four considered dislocation types. The 
lowest Peierls stress value, obtained for the 30° dislocation, is of the order of 1.5 xl0~ 2 MPa. 
This value is ~ 3 orders of magnitude smaller than the shear modulus, which is consistent 
with the typical discrepancy between the ideal shear strength and actual yield stresses in 
crystals [lj]. 
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The above results were obtained using the SWF model, which, while providing a good 
description of the shear modulus for hep solid 4 He, significantly underestimates a recent 
experimental estimate for the SFE, (0.07±0.02) mJ/m 2 |24j. In order to explore the possible 
influence of the SFE value on the dislocation mobility, we repeat the Peierls stress calcula- 
tions for the case in which the equilibrium dissociation widths of the four dislocation types 
is reduced by a factor 10, which is the ratio between the experimental and theoretical SFE 
values. To this end, we apply an additional shear stress component, whose direction in the 
glide plane is perpendicular to the total Burgers vector of the dislocation. This stress com- 
ponent, known as Escaig stress j9j, does not produce a force on the dislocation as a whole, 
but mimics a situation with a different SFE value. Using an Escaig stress of 0.4 MPa, we 
recompute the Peierls stress values for the four dislocations types in the basal plane. The 
results are shown in the fourth row of Table HT1 The effect of an increased effective SFE value 
does not significantly affect the Peierls stress values of the model. This is consistent with 
earlier PN calculations in metals, in which the Peierls stress was not found to be sensitive 
to dissociation width 14|. 



Finally, we examine our results in the context of recent experiments considering the 
macroscopic behavior of solid 4 He at low temperatures. In the observation^, [|| of non- 
classical rotational inertia (NCRI), interpreted as a signature of superfluidity, the crucial 
role of crystal defects and disorder seems firmly established. Specifically, the behavior of 
dislocations has attracted a particular interest after the discovery of an unexpected increase 
of the shear modulus that shows the same temperature and 3 He impurity concentration de- 



pendence as the original 
Granato-Lucke theory 



NCRI observations [6-8]. Inspired by the continuum-elasticity based 



25] , it has been hypothesized that this stiffening is a consequence of 



a change of mobility of a network of dislocations. This network is thought to be pinned by 
3 He impurities at lowest temperatures, while it becomes mobile under warmer conditions. 
Analyzing the dislocation mobility results of our model, it is interesting to observe that our 
lowest Peierls barrier is about 20 times larger than the shear stresses of ~ 700 Pa reached 
in recent experiments This suggests that intrinsic lattice resistance is an essential factor 
when it comes to the mobility of dislocations on the basal plane in hep solid 4 He. Indeed, 
at the stress levels reported in these recent experiments, such dislocations would not be 
expected to be mobile, not even in the absence of any pinning centers. Moreover, it is not 
expected that the Peierls stress varies significantly as a function of temperature below 0.1 K, 
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at which the stiffening is observed, given that finite temperature path-integral Monte Carlo 
calculations of several properties do not show a significant temperature dependence below 1 
K j^. In this context, a satisfactory explanation for the observed elastic stiffening in hep 
solid 4 He must involve intrinsic mobility issues. 

In summary, we employ a hybrid continuum-atomistic approach, based on the Peierls- 
Nabarro model and the shadow wave function formalism, to obtain direct insight into the 
intrinsic structural and mobility properties of dislocations in hep solid 4 He at zero temper- 
ature. In addition to providing key information concerning the elastic properties of this 
prototypical quantum crystal, the results reveal a significant lattice resistance to dislocation 
motion. Analyzing our results in the context of the proposed dislocation-pinning interpre- 
tation of the similarity between the NCRI and elastic stiffening phenomena suggests that 
intrinsic lattice resistance is an essential factor when it comes to the mobility of dislocations. 
The proposed interpretation, which entirely ignores this element, may therefore not provide 
a satisfactory explanation for the observed elastic stiffening in hep solid 4 He. 
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Table I: Elastic constants of the SWF mode/ for hep 4 He within the SWF model, in units of MPa. 





C33 


C44 


C12 


C13 


60.8±0.8 


77.9±0.8 


17.1±0.8 


34.4±0.8 


14.4±0.8 
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Figure 1: Degrees of freedom in the PN model, a) Atomic rows on one side of slip plane (dashed 
line) are labeled by index i. Distance between adjacent rows in defect-free crystal is Ax. b) 
Displacements of rows with respect to those on other side of slip plane are described by vector 5{ . 




Figure 2: Fourier series representation of GSF surface on the basal plane of hep 4 He as modeled by 
the SWF model, a) GSF energy (in mJ/m 2 ) as a function of two-dimensional displacement 5 (in 
units of bp = ^\/3a = 2.0996 A), b) Contour plot of GSF surface reflects symmetry of basal plane 
of hep surface. 
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Figure 3: Optimized displacement profile <5j for screw dislocation on basal plane: edge component 
(open circles) and screw component (filled squares) are measured in units of the Burgers vector 
b = a = 3.63668A. Dashed lines indicate positions of partial dislocations. 

Table II: Core structure and Peierls Stress for four dislocations on the basal plane for the SWF 
model for solid hep 4 He. 





Screw 


30° 


60° 


Edge 


Partial separation (A) 


91.3 


109.1 


119.7 


127.3 


Peierls Stress (MPa) 


0.41±0.01 


0.015±0.01 


0.23±0.01 


0.08±0.01 


Peierls Stress corrected SFE (MPa) 


0.48±0.01 


0.045±0.01 


0.39±0.01 


0.08±0.01 
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